Bose Einstein Condensation of incommensurate solid 4 He 
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It is pointed out that simulation computation of energy performed so far cannot be used to 
decide if the ground state of solid 4 He has the number of lattice sites equal to the number of 
atoms (commensurate state) or if it is different (incommensurate state). The best variational wave 
function, a shadow wave function, gives an incommensurate state but the equilibrium concentration 
of vacancies remains to be determined. In order to investigate the presence of a supersolid phase we 
have computed the one-body density matrix in solid 4 He for the incommensurate state by means of 
the exact Shadow Path Integral Ground State projector method. We find a vacancy induced Bose 
Einstein condensation of about 0.23 atoms per vacancy at a pressure of 54 bar. This means that 
bulk solid 4 He is supersolid at low enough temperature if the exact ground state is incommensurate. 

PACS numbers: 67.80.-s 
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Experiments by Kim and Chan0, Q give evidence for 
non classical rotational inertia of solid 4 He, one hallmark 
of the supersolid state of matter |3j. These results have 
generated large interest because this would be a novel 
state characterized, in a bulk sample, by spontaneous 
broken translational symmetry and by a suitable rigid- 
ity of the phase of the wave function (wf) giving rise 
to superfiuid properties. The standard mechanism for 
this rigidity is the presence of Bose Einstein condensa- 
tion (BEC). Such state with BEC was suggested long 
ago0 Q as a possibility for a quantum solid of boson 
particles. Early theoretical works^ 0, 0, were based 
on simplified models so that it was not possible to draw 
specific predictions for solid 4 He. Powerful simulation 
methods have been applied to study a realistic model 
of solid 4 He in the last two years. Path Integral Monte 
Carlo (PIMC) has been applied to study crystalline He 
at a finite temperature and the authors conclude that the 
superfiuid fraction p s at T = 0.2 K is zero0 and that|| 
there is no off diagonal long range order (ODLRO), i.e. 
there is no BEC, at T = 0.2 K and 0.5 K. On the other 
hand the ground state of crystalline He has been studied 
by variational methods based on shadow wave function 
(SWF) and this study shows that ODLRO is present |2| 
for a range of densities above melting with a rather small 
value of BEC fraction no ~ 0.5 x 10~ 5 at melting density. 

One important point to mention is that in these PIMC 
and SWF computations the crystal is commensurate in 
the sense that the number M of lattice sites is equal to 
the number N of par ticles, i.e. M = N. One finds in 
the liter ature0, IB 113 statements that it is certain that 
the ground state of solid 4 He is commensurate because 
microscopic computations [ill have shown that a va- 
cancy increases the energy of the system by at least 15 K 
and by an even higher value in the case of an interstitial. 
The first purpose of this letter is to present a critical dis- 
cussion of such statements and to examine if the ground 
state of crystalline 4 He is commensurate (i.e. N = M) 
or if it is incommensurate (N ^ M) in the sense that 



the lattice parameter inferred from bulk density measure- 
ment differs from the one deduced from Bragg scattering. 
Present experiments 13] do not give evidence for vacan- 
cies at low T but new measurements seem needed to put 
a stringent bound on ground state vacancies. The nature 
of the ground state, commensurate (C) or incommensu- 
rate (I), is a very important point and a phenomenolog- 
ical theory^] has shown that the low T properties of 
crystalline 4 He would be strongly modified should the 
ground state be incommensurate. Our conclusion will be 
that the microscopic computations of solid 4 He present in 
literature do not allow to infer if the ground state is C or 
I. Since the presence of vacancies in the ground state can- 
not be excluded it is important to study their properties, 
in particular we study if there is BEC induced by vacan- 
cies. This has been studied variationally^5| and here we 
present an exact computation by the projection method 
SPIGS ^2 that confirms a vacancy induced BEC. 

In the first place we notice that the computations [ill 
ll2T ] at the basis of the estimate of the formation energy of 
a vacancy are actually based on the computation of the 
ground state energy of two different systems. To be spe- 
cific let us consider the SPIGS computation. The method 
is based on the application of the imaginary time evolu- 
tion operator exp{— tH} to an assumed trial function, on 
a SWF in the case of SPIGS, and on a splitting of this op- 
erator exp{— tH} = [exp(— -pH)] p which gives rise to a 
path integral of linear polymers. When r is large enough 
a sampling of the exact ground state of the system is ob- 
tained. Notice that at no stage of the computation, either 
in the SWF or in the projection procedure, equilibrium 
sites of the solid are introduced but the crystalline or- 
der, if stable, arises as spontaneous broken symmetry. 
With this method two computations are performed, one 
in which the number N of particles is equal to the num- 
ber M of lattice sites which fit in the simulation box and 
satisfy the periodic boundary conditions (pbc) and one 
in which N = M — 1. In the second case the simula- 
tion shows that the local density continues to have M 
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FIG. 1: Energy per particle as function of the evolution in 
imaginary time r of a SPIGS simulation of fee solid 4 He in 
a box with M = 108. Filled circles: case C (N = M) at 
p = 0.03lA _3 ; open circles: case I (N = M — 1). The largest 
value of r corresponds to 14 projections. Statistical errors are 
below symbol size. 



maxima with essentially the same degree of crystalline 
order, as measured by the height of the Bragg peaks, as 
in the case N = M. This means that in the second case 
the crystalline order is stable with one mobile vacancy 
and such a state is I. The SPIGS computation in the two 
cases gives a converging energy and one example of the 
evolution of E as function of r is shown in Fig.l for the 
fee lattice. Case C is for N = M= 108 at p = 0.03lA~ 3 
and case I is for N = M — 1 = 107. Starting from a 
fully optimized SWF0 just few projections are enough 
to get convergence in both cases. As inter-atomic po- 
tential we have used a standard Aziz potential|l7|. the 
time step S = t/P is (80 K) _1 and the pair-product 
approximation 18] has been used for the imaginary time 
propagator. Both wfs are non negative so both compu- 
tations produce ground state energies but the value Ei 
is slightly larger than Ec- Since we are comparing the 
energy for two different choices of N it makes no sense 
to minimize the energy, both values represent a ground 
state energy of two periodically repeated small systems. 
We conclude that this kind of computation cannot be 
used to determine if the ground state of bulk solid 4 He is 
C or I unless one is able to extrapolate these finite size 
results to the bulk limit also taking into account the ef- 
fects of the pbc. The difference Ej—Eq has been used|p2T] 
to estimate the formation energy of an extra vacancy in 
bulk under the hypothesis of non interacting vacancies 
but this is only a derived quantity [lsj. 

The fact that the constraints imposed by simulation 
of a small periodically repeated system do not allow to 
determine if the ground state of bulk solid 4 He is C or 
I is not surprisin g. In fact also in the classical case it 
is well recognized |20j that no direct computation of the 
equilibrium concentration X v of vacancies of a solid at 
finite temperature is available yet exactly for the same 
reasons as in the quantum case: due to finite size of any 
system that can be simulated and to commensurability 



effect between crystal lattice and simulation box the crys- 
tal cannot achieve its true equilibrium concentration of 
vacancies ll^. X v in the classical case has been obtained 
only indirectly by a statistical thermodynamics analysis 
of an extended system. In a similar way we can expect to 
get light on the nature C or I of the ground state of bulk 
crystalline 4 He only by considering the wf of an extended 
system, not of the one which is simulated. 

In the framework of variational theory the wfs fall in 
two categories. One category is a wf which contains as a 
factor one-body terms which explicitly break the transla- 
tional symmetry. Such wf has great difficulty in describ- 
ing a vacancy and in fact we are not aware of any such 
computation. In addition the nature C or I is built into 
the wf by construction. Lastly such wfs with explicitly 
broken translational symmetry give a worst energy[16| 
than the one given by wfs of the next category. This sec- 
ond category represents translationally invariant wfs for 
which the crystalline state arises as spontaneous broken 
symmetry. One such wf is the time honoured Jastrow 
wf. The other one is SWF which presently gives the best 
representation^^ of the ground state of solid 4 He in the 
sense that it gives the lowest energy. 

An important point is that both these translationally 
invariant wfs give a ground state with a finite concentra- 
tion of vacancies and with BEC. In the case of a Jastrow 
wf this was shown by Chester and we briefly repeat 
here the argument. Consider a Jastrow wf of a very large 
system of TV particles in volume V 
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As noticed long ago 21] computation of averages with 'I'j 
and the normalization have a straightforward inter- 
pretation in classical statistical mechanics: coincides 
with the normalized probability distribution in configu- 
rational space of M classical particles at inverse temper- 
ature /3* = 1/ksT* and interacting with a pair potential 
v*(r) such that (3*v*{r) = u(r). is the canonical con- 
figurational partition function of this classical system and 
its logarithm is proportional to the excess Helmohltz free 
energy. It has been provedj^ that \t j has a finite BEC 
fraction but it is also known that the equivalent classical 
system corresponding to is a crystalline solid, when 
the density is large enough, and this solid has a finite con- 
centration of vacancies. For a classical system the fact 
that a solid in equilibrium at a finite temperature has a fi- 
nite concentration X v = {M —N)/N of vacancies, where 
M is the number of lattice sites, even if a single vacancy 
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has a finite cost of local free energy derives from the gain 
in configurational entropy when the number A4 — Af of 
vacancies is proportional to A/"|13- Another way of ex- 
pressing this is that the configurational partition function 
Q/tf of this equivalent classical system has contributions 
from different pockets in configurational space, from a 
pocket fio in which the positions {r^} of the particles cor- 
respond to vibrations around the equilibrium positions of 
the commensurate Af — AA lattice but also from pockets 
ill, VI2, ... corresponding respectively to AA = Af + 1, 
i.e. a state with one vacancy, to A4 = Af + 2 and so on. 
It turns out [201 that the overwhelming contribution 
to Qjv is associated with pockets Slx-Af with a macro- 
scopic number AA — Af of vacancies. These observations 
have an immediate interpretation in the quantum case: 
the wf P of a bulk system is describing at the same 
time states with no vacancies but also with vacancies 
and the overwhelming contribution to the normalization 
constant Q^f derives from the pockets corresponding to 
a finite concentration of vacancies. The simulation of a 
small system of N particles with pbc is mimicking ex- 
pectation values of the quantum Hamiltonian of the ex- 
tended system in a restricted pocket in configurational 
space, for instance the pocket flo of the commensurate 
state or the pocket fli of the state with one vacancy de- 
pending if N = M or N = M-l. Notice that in a Monte 
Carlo (MC) computation the normalization constant (J2J 
is never computed explicitly but averages are implicitly 
normalized to the set of configurations that are explored 
in the MC simulation, i.e. to the pocket fl or fli that 
one has implicitly chosen at the start of the computation 
by choosing N = M or N = M — 1. If we try to estimate 
the ground state energy per particle eg = Eg/Af of a 
truly macroscopic system the answer is clear as long as 
the concentration of vacancies is small so that they can 
be considered as independent: If eo = Em=n/N is the 
energy per particle from simulation of the C state and 
E\ = Neo + Ae v the total energy from simulation of the 
I state with one vacancy, the inferred ground state energy 
of the extended system is 

e-G = eo + X v Ae v (3) 

where X v is the average concentration of vacancies that 
should be obtained from an analysis of Qjsf of the ex- 
tended system. This is the true variational energy of *$>j 
and not eo- Notice that eo will differ from eo only by a 
very small amount if X v is well below 1%. 

At present the best variational representation of solid 
4 He is given by a SWF0. A SWF ^swf contains not 
only explicit correlations between {r^} like in Q but also 
indirect correlations built via subsidiary variables, one 
for each atom. Also ^ 2 sy/p has a classical interpretation, 
in fact the normalization of ^! 2 g WF coincides with the 
configurational partition function of a classical system of 
suitable flexible triatomic molecules |24|. For this equiva- 
lent classical system the concentration X v of vacancies is 



finite since in the previous argument it makes not differ- 
ence that the "particles" are mono-atomic or are molec- 
ular species. We conclude that also #swf describes a 
quantum solid with vacancies in it and the ground state 
energy of an extended system is given by eal5 ( 27 | . 

Next question is how big is X v . An estimate of X v 
for a Jastrow function has been performed some years 
ago |26| but unfortunately gives an unrealistic repre- 
sentation of solid 4 He because VPj gives a much too large 
localization of atoms. For ^swf X v is not known and 
this is a priority computation for the future. 

We consider now the exact ground state as given by 
SPIGS 12] . The projection maps the quantum prob- 
lem into an equivalent classical problem of flexible lin- 
ear polymers with the number D of monomers equal to 
D = 2P + 1 where P is the number of projections. Such 
classical system has vacancies for any finite P but the 
concentration X v of vacancies might vanish in the limit 
P — > 00. Only a study of X v as function of P will be 
able to say if ground state vacancies are present or not 
in the exact ground state of solid 4 He. 

Given that vacancies might well be part of the ground 
state of solid 4 He we present a microscopic calculation 
of the one-body density matrix pi in presence of va- 
cancies with the exact SPIGS method. Computation of 
pi(r, f) is performed by cutting one of the linear poly- 
mers at the central monomer and sampling the distance 
f — f I of the two cut ends. Notice that in SPIGS com- 
putations, contrary to the case of PIMC, no exchange 
moves between polymers have to be performed and this 
is a great advantage due to ergodicity problems arising 
from exchange moves. In any case calculation of p\ by 
means of SPIGS in the solid is computationally very in- 
tensive due to the low relaxation (to get converged results 
one needs more than 10 7 Monte Carlo steps) and to the 
large number of degrees of freedom (i.e. coordinates of 
monomers) as the imaginary time evolution r becomes 
greater. There is also the necessity to compute pi as a 
function of r in order to control the convergence. We 
have worked with a box which fits M — 108 fee lattice 
sites and with TV = 107 so that we have one vacancy out 
of 107 atoms which corresponds to X v = 0.93%. In Fig[21 
we show the results for three values of r|25| as well as 
the variational SWF results at density 0.031 A -3 , which 
corresponds to a pressure of about 54 bar. p\ has been 
computed with f — r* along the nearest neighbor direc- 
tion ([110] in fee). One can see that p\ develops a plateau 
for distances greater than about 5 A and this is a signa- 
ture of BEC. We have estimated the condensate fraction 
by averaging the plateau for distance greater than 5.5 
A and the values are shown in Fig0 The value of uq 
is only weakly dependent on r and similar to the SWF 
results. With SWF it is known that no for the hep lat- 
tice is very similar to the fee one 15] . pi in the plateau 
region has oscillations and the maxima correspond to 
multiples of the nearest neighbor distance. Since pi is 
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FIG. 2: One-body density matrix pi computed with SPIGS 
in fee solid 4 He at p = 0.03lA -3 for different imaginary time 
evolutions r: (a) r = 0.025 K _1 ; (b) r = 0.075 K _1 ; (c) r = 
0.125 K" 1 ; (d) pi computed with SWF with N = M-l = 107 
(continuous line) and with N — M — 1 = 191 (dashed line). 
Dotted lines represent statistical uncertainty. 



proportional to the probability to destroy a particle at r 
and to recreate it at f 1 one can interpret the process giv- 
ing origin to ODLRO in terms of a sequence of jumps of 
atoms which make use of the vacancy. This process is dis- 
tinct from the vacancy-interstitial pairs that were found 
to be important for the commensurate state in a SWF 
computation!^. The simulated system is small and with 
SWF we have computed pi also for a larger system with 
N = M - 1 = 191 which corresponds to X v = 0.52%. 
One sees from Fig|2 the persistence of the oscillations 
around a finite plateau and it is found that no roughly 
scale with X v . We are confident that this remains true 
also for SPIGS so that we estimate from no and the value 
of X v of our computation that there is a condensate frac- 
tion of about 0.23 4 He atoms per vacancy at the pressure 
of about 54 bar. Therefore vacancies are very efficient in 
inducing BEC. Using, as an order of mag nitude, Tbec of 
an ideal gas with the effective mass 28] m* = 0.35 mu e 
we get T BEC ~ 11.3 x (X„) 2/3 ; for example T BEC ~ 0.2 
K for X v = 2.3 x 10 -3 . Therefore we expect supersolidity 
at low T in bulk solid 4 He if vacancies are present either 
as part of the ground state or as non equilibrium effect. 

In conclusion we have shown, on the basis of an ex- 
act microscopic theory of solid 4 He, the SPIGS projector 
method, the presence, at the same time, of spatial order 
and of BEC when a finite concentration of vacancies is 
present at T — K i.e. if the ground state is incommen- 
surate. Based on the argument by LeggettQ this system 
would show non classical rotational inertia effects. In 
addition we have shown that the question if the ground 
state of bulk solid 4 He is C or I is still undecided but we 
noticed that the ground state is I for the best variational 



wf. It remains an open problem the quantitative evalu- 
ation of the concentration of vacancies X v for the SWF 
and the study of what happens to X v under projection 
with the SPIGS method. 
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